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Abstract Within this decade gravitational wave detection will open a new 

observational window on the Universe. Advanced ground based interferometers 

covering the kHz frequency range will be online by 2016, and it is foreseeable 

CN ■ the announcement of a first detection within five years. At the same time, 

a worldwide effort of detecting low frequency waves (in the nHz regime) by 

/^ ■ timing ultra-precise millisecond pulsars is rapidly growing, possibly leading to 

a positive detection within this decade. The mHz regime, bridging these two 

windows, is the realm of space based interferometers, which might be launched 

rS^ [ in the late 20s. I provide here a short overview of the scientific payouts of 

r^. gravitational wave astronomy, focusing the discussion on the low frequency 

^ I regime (pulsar timing and space based interferometry) . A detailed discussion 

H ■ of advanced ground based interferometer can be found in Patrick Brady's 

t^ ' contribution to this proceedings series. 

C^ ! 

Keywords Relativity • Gravitational waves • Black hole physics • Pulsars: 

general 

^^O ' 1 Introduction 

l> 

^^ I Our current knowledge of the Universe comes primarily from the observa- 

\l ■ tion of photons released in the physical processes driving the formation and 

evolution of stars and galaxies. However, black holes (BHs, among the most 
exciting astrophysical objects) are 'dark' in nature, and can be probed in a 
clean way only by directly capturing dynamical modification of the spacetime 
metric induced by their accelerated motion and radiated away at the speed 
of light. Those are better known as gravitational waves (GWs). GWs are not 



5-H ■ exclusively produced by BHs, but by any accelerated mass with a time varying 
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mass quadrupole moment, however they are so weak, that the only potentially 
observable sources involve astrophysical compact objects. Among the most 
appealing sources are: catastrophic events like supernovae (SNe) explosions 
and gamma ray bursts (GRBs) leaving behind a compact remnant, fast ro- 
tating neutron stars (NSs) with some irregularity in their mass distribution, 
and, in particular, compact binaries. The latter class of sources include any 
combination of stellar compact objects (white dwarfs (WDs), NSs and BHs) 
and systems involving massive and supermassive BHs (to which we will gen- 
erally refer as MBHs in the following): extreme mass ratio inspirals (EMRIs) 
of stellar compact objects into a MBH, and MBH binaries. 



The signals emitted by the aforementioned sources have so far silently 
been lost in space, however, within this decade the detection of GWs may 
indeed be a reality, opening a completely new window on the Universe. While 
signals coming from compact stars and binaries fall in the observational do- 
main of operating and planned ground based interferometers (such as the 
advanced LIGO [J, VIRGO 0, KAGRA 0, and the proposed 'third gener- 
ation' Einstein Telescope (ET) [4]), MBH binaries and EMRIs are expected 
to be among the primary actors on the upcoming low frequency stage, where 
the 10"** — 10~^Hz window is going to be probed by spaceborne interferom- 
eters like the Laser Interferometer Space Antenna (LISA) [S], or the recently 
proposed European development cLISA/NGO [6][3- At even lower frequen- 
cies, precision timing of an array of millisecond pulsars (i.e. a pulsar timing 
array, PTA) provides a unique opportunity to get the very first low-frequency 
detection. The European Pulsar Timing Array (EPTA) [5] , the Parkes Pulsar 
Timing Array (PPTA) [3] and the North American Nanohertz Observatory for 
Gravitational Waves (NANOGrav) [TU], joining together in the International 
Pulsar Timing Array (IPTA) jll| , are constantly improving their sensitivity in 
the frequency range of ~ 10~^ — 10~^ Hz, and the planned Square Kilometer 
Array (SKA) fl2| will provide a major leap in sensitivity. 



In this contribution I will focus on the most relevant GW astrophysical 
sources, meaning that I will not touch the possibility of detecting cosmologi- 
cal backgrounds from e.g., inflation, phase transitions and defects in the early 
Universe and cosmological strings. An overview of the latter can be found in 
[IB] . Among all astrophysical sources particular attention will be devoted to 
binaries of all kind, which are the one we understand better and that will prob- 
ably carry the most valuable and cleanest astrophysical information. Excellent 
extensive reviews on the subject can be found in P^IT5] . After a short general 
introduction about the principles of GW generation and detection (Section [2]), 
I will only briefly touch the ground based interferometer science in Section 3. 
Section 4 is devoted to space-based interferometers and Section 5 to PTAs. 
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2 Gravitational waves: generation and detection principles 

2.1 Gravitational wave generation 

The existence of GWs was one of the first predictions of Einstein's General 
Relativity (GR) [TB]. However, universal acceptance of their existence came 
only in the '80s, thanks to the excellent agreement between theoretical pre- 
dictions [T7] and observations of the Hulse- Taylor binary pulsar [TS] . GWs are 
ripples in the fabric of the spacetime propagating at the speed of light. Ac- 
celerating masses generate GWs pretty much in the same way as accelerating 
charges generate electromagnetic waves. By expanding the mass distribution 
of the source into multipoles, conservation laws enforce GWs coming from the 
mass monopole and mass dipole to be identically zero, so that the first con- 
tribution to GW generation comes from the mass quadrupole moment Q. The 
GW amplitude is therefore proportional to the second time derivative (accel- 
eration) of Q. Moreover, energy conservation enforces the amplitude to decay 
as the inverse of the distance to the source, D. A straightforward dimensional 
analysis shows that the amplitude of a GW is of the order of [TS] 

, GlcfQ 

"=^^^- ^^) 

In order to generate GWs we therefore need accelerating masses with a time 
varying mass-quadrupole moment. The prefactor G /c^ implies that these waves 
are tiny^ so that the only detectable effect is produced by massive compact as- 
trophysical objects, as enumerated in the Introduction. A detailed derivation 
of the GW properties is beyond the scope of this review, I shell just mention 
that GWs are transverse ,i.e. they act in a plane perpendicular to the wave 
propagation, and (at least in GR) have two distinct polarizations, usually re- 
ferred to as /i+ and h^ (see cartoon in [14]). Their quadrupolar nature implies 
that perpendicular directions experience 'opposite' squeezing, so that if the 
GW propagates in the z direction, the metric shrinks in the x direction while 
expanding in the y direction and viceversa in a oscillatory fashion. 

Using equation ^ we can infer approximate GW amplitudes for typical 
sources. Consider, for simplicity, a circular BH binary system characterized by 
a Keplerian frequency Jk = //fl chirp mass M = M^^^M^^^/{Mi + Ms)^/^ 
(Ml > AI2 are the individual masses of the two BHs), and distance to the 
observer D, equation ([1} results in a typical strain of 
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The maximum emission frequency is approximately given by twice the Kep- 
lerian frequency of the binary at the innermost stable circular orbit (ISCO, 



^ A circular binary emits GWs at a frequency which is twice the Keplerian binary fre- 
quency, whereas an eccentric one emits a more complicated wave covering the spectral range 
n/j., where n is an integer index. 
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binary separation a — QGAIijc^ — iRs, where Rs is the Schwarzschild radius 
of the more massive object in the binary) 

/isco ~ lO'Hz^. (3) 

iVi0 

Equations © and ^ are indeed vahd for any binary of compact objects, with 
one important difference: when the binary involves 'extended objects' (hke NSs 
or WDs), the maximum frequency is given by the 'contact' separation of the 
two bodies. For example, the maximum WD- WD binary frequency is roughly 
Q.lHz, beyond which, the two body merge in a single object. 

Plugging in the mass quadrupole moment for a rotating ellipsoid with 
ellipticity e into equation ([T]) gives instead 

-2ef D Y^ f f V f e\ f Q 



where we normalized to typical values of a rapidly spinning NS. On the other 
hand, predicting the waves coming from collapses and explosions is much more 
difficult and requires detailed 3D simulations of the process. A fair bet is to 
assume that the explosions involve the non-axisymmetric acceleration of some 
substantial fraction of Q, to a speed which is a fraction of c. This automatically 
sets the range of amplitudes and frequencies of interest to h ^ 10^^'^ — 10^^'' 
and / ^ 10^ — lO^Hz respectively (see [19] for an excellent review). 



2.2 Gravitational wave detection 

The quantity h given in equation [T] is usually referred to as 'strain', and rep- 
resents the relative stretch of the metric. If two test masses are placed at a 
distance L, a passing GW would result in a oscillatory relative change of their 
distance AL/L ^ h: the measurement of this relative change is at the basis of 
GW detection. The most promising GW detectors are laser interferometers. In 
these devices, the laser beam is split in the two arms of length L, and travels 
forth and back in the arms before being recombined. A passing GW causes a 
AL = L2 — Li ^ hL that can be observed in the interference pattern of the 
recombined light. Earth based interferometers are limited in length (few Km) 
and cannot overcome the noise generated by gradients in the gravitational field 
of the Earth. They can therefore reach sensitivities of the order of ft. ~ 10"^'^ 
only in the kHz frequency range (See, for reference, the Adv. LIGO sensitivity 
curve in Figure[T]). Stellar mass objects are therefore their main targets. To in- 
vestigate lower frequencies, it is necessary to go to space. Here, the possibility 
of sending laser signals between spacecrafts orbiting at large separations from 
each other, together with the absence of gravity gradients, allows to probe the 
mHz frequency regime, with a strain sensitivity of the order of ft ~ 10^^" (Fig- 
ure [Ij. This is the realm of several classes of astrophysical objects including 
wide stellar binaries, EMRIs and MBH binaries. 
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Precision timing of millisecond pulsars (MSPs) provide an alternative (though 
equivalent) way to detect GWs (see, e.g., pUll^l^ l. MSPs are the most sta- 
ble natural clock in the Universe [23], and the time of arrival of their pulses 
can be currently predicted with a precision of the order of ^ 100ns (see, e.g., 
[24)). If the pulsar and the Earth are in 'free-fall' in a spacetime perturbed by 
the passage of a GW, this will leave a characteristic fingerprint in the time of 
arrival of the radio pulses. As a matter of fact, the passing wave is 'changing 
the distance' between the pulsar and the Earth, which are now the 'two test 
masses' of our detector. The signal can be recovered by correlating the time of 
arrivals of an array of MSPs, i.e. a PTA, which can be (by analogy) thought 
as a 'galactic scale interferometer'. In this latter case, the sensitivity window 
is set by the duration of the timing experiment (typically several years) and 
by the sampling rate of the pulsars (typically once every few weeks) , resulting 
in a range ~ 10"^ < / < 10'^ Hz. 



3 Ground based interferometers 

The status of ground based interferometers and the expected Science deliv- 
ered by their upcoming second generation is presented by Patrick Brady in his 
contribution to this proceeding series, here I just provide a general overview. 
Given their sensitivity curves (Figure [T|), advanced ground-based interferom- 
eters are sensitive to stellar mass sources. Their main goal is the detection of 
inspiralling compact objects (NS-NS, NS-BH, BH-BH), and I report in table 
(P) the expected detection rates, adapted from [27]. The large uncertainties in 
the figures are due to the poor observational constrains (we know only a hand- 
ful of NS-NS systems which are going to coalesce in an Hubble time, [28 ), and 
to the many unknowns affecting theoretical population synthesis models (e.g., 
the physics of common envelope evolution f5^). GW observations will likely 
help to better understand binary evolution, possibly constraining the detailed 
physics of mass transfer and common envelope physics [30] ) . NS deformation 
during the late inspiral of a NS-NS or NS-BH binary will inform us about the 
stiffness of the NS equation of state, probing the behavior of matter in extreme 
degenerate conditions (see, e.g., ^). Coalescing NS-NS and NS-BH are also 
believed to be the engines powering short GRBs (see, e.g., [351), and coincident 
multimessenger (GW plus electromagnetic) detections will therefore provide 
the ultimate test for such hypothesis. Rotating NSs in the Galaxy are another 
potentially interesting source [5S]. Even though the strain is fairly small (see 
equation (|3|)), many wave cycles will accumulate in a prolonged observation, 
resulting in an effective amplitude of the order of hy/N, where N — T/P, 
being P the period of the pulsar and T the observation time (order of years). 
As shown in Figure [1] Adv. LIGO might detect Galactic NSs spinning at kHz 
frequencies with e w 10~^. Collapsing and exploding objects (SNe, GRBs) 
are hard to model (See [19j for a review), and different proposed scenarios 
result peculiar GW signatures roughly falling in the cyan area in Figure [1] 
Detection of GW bursts from SN explosions will therefore help to understand 
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Table 1 Adv. LIGO detection rates, adapted from |27l . From the left to the right numbers 
arc for pessimistic, realistic and optimistic rates (see | 27| for details). The large range stems 
from the uncertainties in the complex physics involved in population synthesis models. 



Source I TViowIyr"^] A^rolyr-l] A^highbr"^] 



NS-NS 
NS-BH 
BH-BH 



0.4 40 400 

0.2 10 300 

0.4 20 1000 



the physical mechanisni driving the collapse and the subsequent catastrophic 
bounce. 



4 Space based interferometers 

Space based interferometers like the proposed eLISA (which will be our default 
design in the following discussion) are sensitive in the mHz frequency band 
(Figure [ij, where a rich variety of sources is expected. The enormous scientific 
payouts related to space based interferometers are far too many to be covered 
in a short review. The reader can only get a sense of it from the following 
discussion, a complete overview can be found in [6j[7]. 

Space based interferometers are the only detectors for which we have 'ver- 
ification sources' |33| . These are known ultra-compact binaries emitting in 
the mHz range with high signal-to-noise ratio (SNR), and will serve as de- 
tector calibrators. Assuming the eLISA design, eight such binaries are known 
to date. Moreover eLISA will be sensitive to the whole population of wide 
WD- WD systems (order 10* sources) populating the Galaxy. Most of those 
adds-up to the overall confusion foreground |M1I55] , but about 3000 of them 
will be individually resolved [Sl[7]. Many resolvable binaries have periods be- 
tween 5 and 10 minutes implying that they must have experienced at least one 
common-envelope phase; their statistics will therefore provide critical tests of 
physical models of common-envelope binary evolution. Although rarer, eLISA 
will also detect a substantial population of short-period NS and BH binaries, 
determining their local merger rate. 

One of the most promising science goals of eLISA is the direct detection of 
MBHs, which appear to be a key component of galaxies. MBHs are ubiquitous 
in nearby massive galaxies (e.g., [36j ) and share a common evolutionary path 
with their hosts, undergoing repeated episodes of merger and accretion |37[ 
155] . Hierarchical models of MBH assembly predict that 10— 100 MBH mergers 
per year occur in the Universe (see, e.g., [5^l26j ). the bulk of which involves 
systems with lO^'M© - lO'^Mg MBHs at 3 < z < 10, falling in the sweetspot 
of the eLISA sensitivity curve. In fact, eLISA will be able to detect the GWs 
emitted by MBH binaries with total mass as small as 10"*^ M© and up to 10* Mq , 
out to a redshift as remote as 2: ~ 20 with a SNR > 10. Redshifted masses will 
be measured to an unprecedented accuracy, up to the 0.1-1% level, whereas 
absolute errors in the spin determination are expected to be in the range 0.01- 
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Fig. 1 Summary of the GW landscape discussed in the text. Black solid curves are detector 
sensitivities as labeled in the figure. For the PTA sensitivity curve we considered an SKA- 
like scenario in which 500 MSPs are monitored for 10 yrs at 50ns precision. In the Adv. In 
the LIGO band, the blue curves represent inspirals of compact stellar mass binaries |14l . 
the dashed magenta line marks the expected strain for rapidly spinning NSs 1251 . and the 
shaded cyan box identifies the area where burst signals from collapsing SNe are expected |19l . 
In the eLISA/NGO band, the red curves represent merging MBH binaries at cosmological 
distances, the dashed-magenta line is a typical eccentric EMRI at z = 0.1 (only the first 
5 harmonics are shown) and the jagged yellow curve is the overall signal from a galactic 
population of WD-WD binaries. In the PTA band, the total signal (green solid line) is 
given by an incoherent superposition of individual MBH binaries (green dots); the brightest 
sources can be individually resolved (blue lines and triangles). The lower-left inset shows 
the inferred cosmological merger rate of MBH binaries as a function of mass for two specific 
models |26[ . eLISA/NGO will probe the low mass part of the distribution (red shade) by 
detecting individual merging systems, whereas PTA will probe the high mass end (green 
shade) of the mass function by collecting signals from very massive systems in their inspiral 
phase. 



0.1, allowing us to reconstruct their cosmic evolution [33]. eLISA observations 
hence have the potential of constraining the astrophysics of MBHs along their 
entire cosmic history, in a mass and redshift range inaccessible to conventional 
electromagnetic observations. 

On smaller scales, eLISA will also offer the deepest view of galactic nuclei, 
exploring the dynamics of stars in the space-time of a Kerr BH, by capturing 
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the GWs emitted by stellar BHs orbiting the central MBH (EMRIs). A 10 M© 
BH inspiralling into a 10'' Mq - 10'^ M© MBH can be detected up to z ~ 0.5 - 
0.7 [6j[7], allowing to explore a volume of several tens of Gpc"^ and to discover 
MBHs in dwarf galaxies that are still elusive to electromagnetic observations. 
Detailed measurement of the MBH mass and spin (up to one part in a thousand 
[ID] ) will enable us to characterize the population of nuclear MBHs in an 
interval of masses where electromagnetic observations are poor, incomplete or 
even missing 1411 . providing information also on their spins. eLISA will also 
measure with equivalent precision the mass of the stellar BH in the EMRI 
event informing us about the mass spectrum of stellar BHs, which is basically 
unknown. 

High SNR measurement of the GW waveforms will allow to test GR in the 
highly dynamical strong field regime of MBHs, including no-hair theorem tests, 
constrains on a putative graviton mass and on alternative theory of gravity in 
general. More details about testing gravity with GW observations are given 
by Emanuele Berti in his contribution to this proceedings series |42) . 



5 Pulsar timing arrays 

By exploiting the characteristic fingerprint left by passing GWs in the time 
of arrival of the radio pulses propagating from the pulsar to the receiver on 
Earth [20ll21l[22] . PTAs are sensitive to the collective signal coming from the 
inspiralling population of supermassive low redshift binaries [^Hi] in the nHz 
frequency band. The overall expected characteristic strain he can be written 
as 

where (PN/dzdAidlnfr, is the comoving number of binaries emitting in a 
given logarithmic frequency interval with chirp mass and redshift in the range 
[A4,A4+dA4] and [z, z+dz], respectively; and h{fr) is the sky and polarization 
averaged version of the strain amplitude given in equation ([1]) , and is given by 

m- 

101/2 i?i(z)^'- ' ^^^ 

where Dl{z) is the luminosity distance to the source. The GW spectrum has 
a characteristic shape he — A{f /lyi^^)^^^^ (green solid line in Figure [1]), 
where A is the signal normalization at / = lyr~i, which depends on the 
details of the MBH binary population only. Therefore, GW detection will 
provide a strong test for the effectiveness of the merger process of galaxies 
and their MBHs at low redshift, and will help constraining the high mass 
end of the MBH mass function. Recent works [46u47y48| set a plausible range 
3 X 10^1^ < A < 4 X 10~i^, the upper limit being already in tension with 
current PTA measurement [¥^[50] . With few more years of observation, PTAs 
have a concrete chance of detecting this signal. Several millions of sources con- 
tribute to it, however, the bulk of the strain comes from few hundred sources 
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only (see, e.g., [5T]) generating a signal at a Ins level or higher, considered 
the ultimate goal for the SKA. Therefore, the signal is far from being a Gaus- 
sian isotropic background |46j ; a handful of sources dominate the strain bud- 
get, and might be individually identified. Interestingly, such systems are far 
from coalescence, and they can still retain much of their original eccentricity 
against GW circularization [5^[551[5i] . Eccentricity measurement of individu- 
ally resolved sources may help in constraining the evolution of MBH binaries, 
testing our current models of their dynamical evolution in star and/or gas 
dominated environments. For bright enough sources (SNRss 10) sky location 
within few tens to few deg^ is possible [5^F5^I?7] (and even sub deg^ determi- 
nation, under some specific conditions [55]). Even though this is a large chunk 
of the sky, these systems are extremely massive and at relatively low redshift 
{z < 0.5), making any putative electromagnetic signature of their presence 
(emission periodicity related to the binary orbital period, peculiar emission 
spectra, peculiar Ka line profiles, etc.) detectable [551150] . 



6 Conclusions 

Gravitational wave astronomy has never been so active as in the last few 
years. The prospect of the first detection from Advanced ground based inter- 
ferometers has triggered a wave of excitement in the astrophysical community, 
pushing the concept of multimessenger astronomy and the development of 
pipelines that combine GW and electromagnetic interactions to unveil the na- 
ture of short gamma ray bursts and supernovae. A worldwide effort of detecting 
low frequency waves (in the nHz regime) by timing ultra-precise millisecond 
pulsars is running in parallel. The enormous advanced in the quality of the 
receivers and in the data reduction techniques is making the goal of achiev- 
ing sub-lOOns timing precision on a bunch of pulsars possible. Current array 
sensitivities are a factor of a few away from theoretical prediction of the nHz 
signal expected from cosmological supermassive black holes, making detection 
within this decade plausible. Looking a step ahead in the future, space based 
interferometer such as LISA or eLISA/NGO will cover the mHz regime, pro- 
viding an almost full coverage of the gravitational wave spectrum from the 
nHz to the kHz. The silent side of the Universe will soon get loud. 



References 

1. Harry G. M., and the LIGO Scientific Collaboration, 2010, Classical and Quantum Grav- 
ity, 27, 4006 

2. Acernese F., et al., 2006, Classical and Quantum Gravity, 23, 63 

3. Somiya K., et al., 2012, Classical and Quantum Gravity, 29, 4007 

4. Punturo M., et al., 2010, Classical and Quantum Gravity, 27, 4002 

5. Bender P.-L. et al., 1998, LISA Pre-Phase A Report, Second Edition, MPQ 233 

6. Amaro-Seoane P., et al., 2012, Arxiv e-prints 1201.3621 

7. Amaro-Seoane P., et al., 2012, Classical and Quantum Gravity, 29, 4016 

8. Fcrdman R. D., et al., 2010, Classical and Quantum Gravity, 27, 084014 



10 Alberto Scsana 



9. Manchester R. N., et al., 2012, ArXiv e-prints 1210.6130 

10. Jenet F. A., et al., 2009, ArXiv e-prints 0909.1058 

11. Hobbs G., et al., 2010, Classical and Quantum Gravity, 27, 084013 

12. Aharonian F., et al., 2013, Arxiv e-prints 1301.4124 

13. Binetruy P., Bohe A., Caprini C. & Dufaux J.-F., 2012, JCAP, 6, 27 

14. Thorne, K. S., Particle and Nuclear Astrophysics and Cosmology in the Next Millenium, 
Proceedings of the 1994 Snowmass Summer Study held 29 June - 14 July, 1994. Edited 
by E.W. Kolb and R.D. Peccei. Singapore: World Scientific, 1995., p. 160 

15. Hughes S. A., 2003, Annals of Physics, Volume 303, Issue 1, p. 142-178 

16. A. Einstein, 1916, Koeniglich Preu§ische Akademie der Wissenschaftcn Berlin, Sitzungs- 
berichte, 688. 

17. Damour T., 1983, Phys. Rev. Lett., 51, 1019. 

18. Taylor J. H. & Weisberg J. M., 1989, Astrophys. J., 345, 434 

19. Ott C. D., 2009, Classical and Quantum Gravity, 26, 3001 

20. Sazhin M. V., 1978, Soviet Astron., 22, 36 

21. HeUings R. W. & Downs G. S., 1983, Astrophys. J., 265, 39 

22. Jenet F. A., Hobbs G. B., Lee K. J. & Manchester R. N., 2005, Astrophys. J., 625, 123 

23. Hobbs G., et al., 2011, MNRAS, 427, 2780 

24. Verbiest J. P. W., et al., 2009, MNRAS, 400, 951 

25. Owen B. J., Lindblom L., Cutler C, Schutz B. F., Vecchio A. & Andersson N., 1998, 
Phys. Rev. D, 58, 4020 

26. Sesana A., Gair J., Berti E. & Volonteri M., 2011, Phys. Rev. D, 83, 4036 

27. Abadie J., et al., 2010, Classical and Quantum Gravity, 27, 3001 

28. Kalogera V., et al., 2004, Astrophys. J. Lett., 601, 179 

29. Dominik M., et al., 2012. Astrophys. J., 759, 52 

30. Gerosa D., Kesden M., Berti E., O'Shaughnessy R. & Sperhake U., 2013, ArXiv c-prints 
1302.4442 

31. Read J., et al., 2009, Phys. Rev. D, 79, 4033 

32. Kochanek C. S. & Piran T., 1993, Astrophys. J. Lett., 417, 17 

33. Stroeer A. & Vecchio A., 2006, Classical and Quantum Gravity, 23, 809 

34. Nelemans G., Yungelson L. R., Portegies Zwart S. F., 2001, A&A, 375, 890 

35. Nissanke S., VaUisneri M., Nelemans G. & Prince T. A., 2012, Astrophys. J., 758, 131 

36. Magorrian J. et al., 1998, AJ, 115, 2285 

37. Volonteri M., Haardt F. & Madau P., 2003, Astrophys. J., 582, 599 

38. Malbon, R. K., Baugh, C. M., Frenk, C. S., & Lacey C. G. 2007, MNRAS, 382, 1394 

39. Sesana A., Volonteri M. & Haardt F., 2007, MNRAS, 377, 1711 

40. Barack L. & Cutler C, 2004, Phys. Rev. D., 69, 2005 

41. Gair J. R., Tang C. & Volonteri M., 2010, Phys. Rev. D, 81, 4014 

42. Berti E., 2013, ArXiv e-prints 1302.5702 

43. Sesana A., Vecchio A. & Colacino C. N., 2008, MNRAS, 390, 192 

44. Sesana A., Vecchio A. & Volonteri M., 2009, MNRAS, 394, 2255 

45. Thorne K. S., 1987, Gravitational radiation., pp 330458 

46. Ravi v., et al., 2012, Astrophys. J., 761, 84 

47. McWilliams S. T., Ostriker J. P., Pretorius F., 2012, ArXiv e-prints 1211.4590 

48. Sesana A., 2012, MNRAS, in press, ArXiv e-prints 1211.5375 

49. van Haasteren R., et al., 2011, MNRAS, 414, 3117 

50. Demorest P. B., et al., 2013, Astrophys. J., 762, 94 

51. Kocsis B. & Sesana A., 2011, MNRAS, 411, 1467 

52. Sesana A., 2010, ApJ, 719, 851 

53. Preto M., Berentzen L, Berczik P. & Spurzem R., 2011, Astrophys. J., 732, 26 

54. Roedig C, Dotti M., Sesana A., Cuadra J. & Colpi M., 2011, MNRAS, 415, 3033 

55. Sesana A. & Vecchio A., 2010, Phys. Rev. D, 81, 4008 

56. Ellis J. A., Siemens X. & Creighton J. D. E., 2012, Astrophys. J., 756, 175 

57. Petiteau A., Babak S., Sesana A. & de Araujo M., 2013, Phys. Rev. D, in press, ArXiv 
e-prints 1210.2396 

58. Lee K. J., et al., 2011, MNRAS, 414, 3251 

59. Sesana A., Roedig C, Reynolds M.T. & Dotti M., 2012, MNRAS, 420, 860 

60. Tanaka T., Menou K. & Haiman Z., 2012, MNRAS, 420, 705 



